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In this proceeding, we discuss recent results on exclusive radiative B meson decays from the Belle Collaboration. 
These decays are sensitive to right-handed currents from Nevf Physics. In particular, we measure time-dependent 
CP violation parameters in ij" — > K^tt^^ and — > K^p^^ decays, using high-statistics data samples collected 
at the T(45) resonance with the Belle detector at the KEKB asymmetric-energy e+e~ collider. With the 
present statistics, these measurements are consistent with the standard model predictions. We also search for 
the radiative decay B'^ — » (pK^-y and report the first observation with a significance of 5.4 tr, including systematic 
uncertainties. 



1. Introduction 

Rare radiative decays of B mesons play an impor- 
tant role in the search for physics beyond the standard 
model (SM) of electroweak interactions. These flavor 
changing neutral current decays are forbidden at tree 
level in the SM, but allowed through the electroweak 
penguin processes as in Fig. 1. Hence, they are sen- 
sitive to non-SM particles mediating the loop (for ex- 
ample, charged Higgs or SUSY particles), which could 
affect either the branching fraction or CP violation. 




Figure 1: Feynman diagram for the radiative b —> s-y de- 
cays, showing the SM loop process with the t-quark con- 
tribution. 

In the SM, the photon emitted from a (-6°) me- 
son in the radiative 6 — s- S7 decays, is predominantly 
right-handed (left-handed). Therefore, the polariza- 
tion of the photon carries information on the origi- 
nal b flavor and the time-dependent CP asymmetry is 
suppressed by the quark mass ratio {2ms /mb) [!]• In 
several models beyond the SM, the photon acquires 
an appreciable right-handed component due to the 
exchange of a virtual heavy fermion in the loop pro- 
cess, resulting in large values of mixing-induced CP 
asymmetries. The same argument holds in any multi- 
body final states B° P^Q^-f, where P° and Q° 
are charge-conjugate states [2] (e.g., i?" Kgp^^, 
(j)Kg'j) [3]. A non-zero value of CP asymmetry 
will be a clear hint of new physics. 



2. Experimental Apparatus 

The Belle detector is a large-solid-angle magnetic 
spectrometer that consists of a silicon vertex detector 



(SVD), a 50-layer central drift chamber (CDC), an ar- 
ray of aerogel threshold Cherenkov counters (ACC), 
a barrel-like arrangement of time-of-flight scintillation 
counters (TOF), and an electromagnetic calorimeter 
(ECL) comprised of CsI(Tl) crystals located inside a 
superconducting solenoid coil that provides a 1.5 T 
magnetic field. An iron flux-return located outside 
the coil is instrumented to detect mesons and to 
identify muons (KLM). The detector is described in 
detail elsewhere [4] . Two different inner detector con- 
figurations were used. For the first sample of 152 x 10^ 
BB pairs, a 2.0 cm radius beampipe and a 3-layer sil- 
icon vertex detector (SVDl) were used; for the latter 
samples, a 1.5 cm radius beampipe, a 4-layer silicon 
detector (SVD2), and a small-cell inner drift chamber 
were used. 



3. Analysis Technique 

At the KEKB asymmetric-energy e+e^ (3.5 on 
8.0 GeV) collider [5], the T(4S') is produced with a 
Lorentz boost of — 0.425 nearly along the z axis, 
which is defined as opposite to the e"*" beam direction. 
In the decay chain T(4S') B^B° /rcc/tag, where 
one of the B mesons decays at time tree to a final 
state /rec, which is our signal mode, and the other de- 
cays at time ttag to a final state /tag that distinguishes 
between and B", the decay rate has a time depen- 
dence given by 



r{At) 



-\At\/T^a 

Atbo 



S sm{AmdAt) 



Acoa{AmdAt) }. (1) 



Here S and A are the CP- violation parameters, Tgo 
is the neutral B lifetime, Arrid is the mass differ- 
ence between the two neutral B mass eigenstates. 
At — tree — ^tag, and the 6-flavor charge q equals -1-1 
(—1) when the tagging B meson is identified as B'^ 
(i?"). Since the and B'^ are approximately at rest 
in the T(4S') center-of-mass system (cms). At can be 
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determined from the displacement in z between the 
/rcc and /tag decay vertices: Ai ~ Az/{Pjc), where c 
is the speed of hght. 



3.1. Flavor Tagging and Vertex 
Reconstruction 

The b flavor of the accompanying B meson is iden- 
tified by a tagging algorithm [6] that categorizes 
charged leptons, kaons, and A baryons found in the 
event. The algorithm returns two parameters: the 
^-flavor charge q, and r, which measures the tag qual- 
ity and varies from r = for no flavor discrimina- 
tion to r = 1 for unambiguous flavor assignment. If 
r < 0.1, the accompanying B meson provides negli- 
gible tagging information and we set the wrong tag 
probability to 0.5. Events with r > 0.1 are divided 
into six r intervals. The wrong tag fractions for the 
six r intervals, wi {I — 1, 6) and possible differences in 
wi between B^ and decays (Awi) are determined 
using high-statistics control samples of semi-leptonic 
and hadronic b ^ c decays [7, 8]. 

The vertex position of the signal-side decay is re- 
constructed from the Kg trajectory in Kgir'^j 
and from the daughters of the in i?° Kgp'^j 
mode, with a constraint on the interaction point. The 
tracks are required to have enough hits in the SVD for 
vertexing. The tag-side B vertex is determined from 
well reconstructed tracks that are not assigned to the 
signal side. 

We determine S and A by performing an unbinned 
maximum-likelihood (UML) fit to the observed At dis- 
tribution. The likelihood fimction is 



CiS,A) = ]Jr^iS,A; At,), 



(2) 



where the product includes all events in the fit. The 
probability density function (PDF) is given by 



V. = (l-ZoOl ^/,7',(At')i?,(Ai,-A0 

j 

+ /olPol(Ai,)- 



d{At') 
(3) 



where j runs over the signal and all background com- 
ponents. The fractions of each component (fj) depend 
on the r region and are calculated on an event-by- 
event basis as a function of the fitted variable. Rj 
is the At resolution function and Pol (At) is a broad 
Gaussian function that represents an outlier compo- 
nent with a small fraction foi- The only free param- 
eters in the final fit are S and A, which are deter- 
mined by maximizing the likelihood function given by 
Eq. 2. We define the raw asymmetry in each At bin 
by {N+ ~ N^)/{N+ + N_), where N+ (iV_) is the 
number of observed candidates with q = +1 (—1). 



4. Time-dependent Analysis of 

This analysis is done in Belle using 535 x 10^ BB 
pairs [9]. Since the time-dependent CP asymmetry is 
not expected to change significantly as a function of 
Kstt'^ invariant mass, we perform two measurements: 
one for B° K*°{^ K^Tr°)j by requiring Mj^o^o to 

lie in the range 0.8 < Mj^o^o < 1.0 GeV/c^, and the 
other for the full range of M^o^o below 1.8 GeV/c^. 

The primary signature of these type of decays is 
the high energy prompt photon. They are selected 
from isolated ECL clusters, with center-of-mass (cms) 
energy in the range 1.4 to 3.4 GeV. The polar angle 
of the photon direction in the laboratory frame is re- 
stricted to the barrel region of the ECL (33° < 0^ < 
128°) for SVDl data, but is extended to the end-cap 
regions (17° < 9^ < 150°) for SVD2 data due to the 
reduced material in front of the ECL. The selected 
photon candidates are required to be consistent with 
isolated electromagnetic showers, i.e., 95% of the en- 
ergy in an array of 5 x 5 CsI(Tl) crystals should be 
concentrated in an array of 3 x 3 crystals and should 
have no charged tracks associated with it. We also 
remove photons from Tr^{r]) — > 77 using a likelihood 
function described in Ref. [10]. Neutral kaons {Kg) 
are reconstructed from two oppositely charged pions 
that have an invariant mass within ±6 MeV/c^ of the 
Kg mass. Neutral pions (tt*') are formed from two 
photons with an invariant mass within ±16 MeV/c^ 
of the 7r° mass. 
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Figure 2: (a) AE distribution within the i\fbc signal slice 
and (b) Mbc distribution within the AE signal slice for the 
whole Mf^o ^0 region. Points with error bars are data. The 
solid curves show the fit results. The dotted curves show 
the qq background contributions, while the dashed curves 
show the sum of qq and BB background contributions. 



The B candidates are identified using two kine- 
matic variables: the energy difference AE = E'^^ — 
-^boam ^^'^ beam-energy-constrained mass Mbc = 



Vi^h^im)^ - (pW)^^ where E^f^^ is the beam en- 
ergy in the cms, and Eg'^^ and p'j^^ are the cms energy 
and momentum, respectively, of the reconstructed B 
candidate. The signal region in AE and Mbc, which 
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is used for the measurements of CP- violating param- 
eters, is defined as -0.2 GeV < AE < 0.1 GeV and 
5.27 GeV/c2 < Mbc < 5.29 GeV/c^. 

After all selections are applied, we obtain 4078 can- 
didates in the AiJ-Afbc fit region, of which 406 are 
in the signal box. The signal yield is obtained from 
an UML fit to the Ai?-Afbc distribution as shown in 
Fig. 2. Figure 3 shows the At distributions of the 
events with 0.5 < r < 1.0 for q = +1 and q = —1 and 
the raw asymmetry. The CP violation parameters for 
KsTT^j for the full K^tt^ invariant mass region as well 
as for the mass region around A'* (892)° are summa- 
rized in Table I. 
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Figure 3; (Top) Proper time distributions for 
Ksn°'y for q = +1 (left) and g = -1 (right) with 
0.5 < r < 1.0. The solid curve shows the total and dashed 
curve shows the signal component. (Bottom) Asymmetry 
in each At bin with 0.5 < r < 1.0. The solid curve shows 
the resuh of the UML fit. 



5. Time-dependent Analysis of 

The first measurement of time-dependent CP asym- 
metry in Kgp^'y mode was performed by Belle 
using 657 x 10^ BB pairs [11]. The advantage of this 
mode is that the B'^ decay vertex can be reconstructed 
from two charged pions from the decays, thus avoid- 
ing the complications and efficiency loss from Kg ver- 
texing. The expected S has opposite sign to that of 

The signal is reconstructed in the decay B^ 
Kgp'^j with p'^ TT^TT" and Kg tt+tt^. The selec- 
tion criteria for high energy prompt photon and neu- 



tral kaons are same as those described in section 4, 
except the photons are required to lie in the bar- 
rel region of the ECL and the Kg invariant mass 
should be within ±15 MeV/c^ of its nominal mass. 



We also reconstruct the 



K^ 



7 decay to 



study the i^TTTr system and to serve as a control sam- 
ple. The K^TT^TT'^ and Kgir^n^ invariant masses are 
required to be less than 1.8 GeV/c^. The B candi- 
dates are selected from the Kgir'^ir'^j sample by re- 
quiring the TT+TT^ invariant mass to lie in the p° re- 



gion, 0.6GeV/c2 < 



< 0.9GeV/c2. Since the 



p is wide, other modes that are not self-conjugate, 
such as K*^ir^j may also contribute. Therefore, we 
first measure the effective CP-violating parameters, 
iScft and ^off , using the final sample and then convert 
them to the CP-violating parameters of P° — > Kgp'^j 
using a dilution factor V, described in Ref. [11]. 

The signal is extracted from an UML fit to the 
A/bc distribution as shown in Fig. 4. The requirement 
-0.1 GeV < A£; < 0.08 GeV is applied. We obtain 
299 events in the signal Afbc region after vertexing. 
Out of these we find a signal of 212 ± 17 events with a 
fraction of 6.0% self-cross-feed (SCF), 53.4 ± 2.6 con- 
tinuum, along with 7.8 K*j, 21.5 other ^^7, and 9.0 
BB events. 
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Figure 4: A-Ihc distributions for B^' — > Kg'jT'^Tv~'y events. 
Points with error bars are data. The curves show the re- 
sults from the r dependent Mbc fit. The dashed and dash- 
dotted curves are the qq and all background. The thin 
curve is the total signal including SCF and the thick curve 
is the total PDF. 



We obtain Scs = 0.09 ± 0.27(stat.)t5] [!^(syst.) and 
Acs = 0.05±0.18(stat.)±0.06(syst.) from an UML fit 
to the observed At distribution. The parameter ScS 
is related to S for Kgp^'j with a dilution factor V = 
Scff/Sj^OpO^, that depends on the K*^tt^ components 
and allows for interference: 

^ ^ imi' + 2Re{FXFB) + F^{K)FBiK)] 

f[\FA\^ + 2ReiFXFB) + \FB\^] ' ^ ' 

where Fa , Pb are photon- helicity averaged amplitudes 
for PO ^ A:gp0(7r+7r-)7 and P" ^ K*^{K^c;T:'^)tt^j, 
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respectively. The factors Fb{K), Fb{K) distinguish 
between K* 



7r+7 and K 

r.0 



from the charged Decay Mode 



7. The phase space 
integral is over the p'^ region. We measure the di- 
lution factor V to be O.SSlo Jg 
mode —>■ K^n^ n^"f using a combination of var- 
ious kaonic resonances with spin > 1 to model the 
Kirn system. By combining it with Scs, we obtain 
SkspOj = 0.11 ± 0.33(stat.)t[J°^(syst.). The fits to 
the observed At distributions and the raw asymmetry 
are shown in the Fig. 5. The fit results are summa- 
rized in Table I. Figure 6 shows the distributions for 
the mTTTT spectrum. 
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Figure 5: Fit projections on the At distributions with (a) 
q = -1-1 and (b) q = —1 for events with r > 0.5. The solid 
curves are the fit while the dashed curves show the back- 
ground contributions. The raw asymmetry as a function 
of At is shown in (c) with a fit curve superimposed. 



200 
180 
160 
140 
120 
100 
80 
60 
40 
20 

3!" 
30 
25 
20 
15 
10 
5 










(a) 


l^-i - r r : ', 1 , , , 1 ' 






(b) 







0.4 0.6 0.8 1 1.2 

in„, (GeV/c^) 



Figure 6: rrinn distributions for (a) K^-k^tx^^ and 

(b) — > K'g-K'^'K~'y . The curves follow the convention 
in Fig. 4. The thin dashed curve is the correctly recon- 
structed B ^ J-i'i (1270)7 signal. 



6. First Observation of ^K^l 
decay 

This decay mode has advantages similar to 
Kgp^^ in the search for new physics. Here the i?" de- 
cay vertex can be reconstructed from the two charged 



Table I Results of the fits to the At distributions. The first 
errors are statistical and the second errors are systematic. 



S 



A 



B° -* K*{m2f-{ 
B° Ks-!v°'y 



-0.32+^;^° ±0.05 



-0.20 ±0.24 ±0.05 
-0.10 ± 0.31 ± 0.07 -0.20 ± 0.20 ± 0.06 
0.11±0.33l2 S^ 0.05 ±0.18 ±0.06 



kaons from the (j) decay. The branching fractions for 
the charged (j)Kj mode and an upper limit on the neu- 
tral mode have already been reported by the Belle [12] 
using 96 X 10^ BB pairs. We search for the neutral 
mode using the full data sample, nearly eight times 
larger than was used in our previous measurement and 
report the first observation with a significance of 5.4 a. 

The signal is reconstructed in the decay £?+ 
(j)K+-/ and B° (jjK^J, with cj) K+R- and 
TT+TT^. The invariant mass of the (j) candi- 
dates is required to be within —0.01 < M^+k- ~ 
< +0.01 GeV/c^, where ni^ denotes the world- 
average (j) mass [13]. The selection criteria for the 
high energy prompt photon and neutral kaons are 
the same as those described in section 4, except the 
photons are required to lie in the barrel region and 
the invariant mass of the tt'^tt" combinations has to 



be in the range 0.482 GeV/c^ < A'U+^- < 0.514 
GeV/c^. The B candidates are selected with a re- 
quirement 5.2 GeV/c^ < Mbc < 5.3 GeV/c^ and 
-0.3 GeV < AE < 0.3 GeV. We define the signal 
region as 5.27 GeV/c^ < A4c < 5.29 GcV/c^ and 
-0.08 GeV < A£; < 0.05 GeV. 

The dominant background is from the continuum 
process, which is suppressed by a requirement on like- 
lihood ratio from event shape variables and the B 
fiight direction. In the B^ mode, some 

backgrounds from b c decays, like D^tt", D'^rj and 
p^ peak in the Afbc distribution. We remove these 
backgrounds by applying a veto on the 4'K'g invariant 
mass. The non-resonant background B K~ Kj, 

which peaks in the Ai?-Mbc signal region, is estimated 
using the 4> mass sideband in data. 

The signal yield is obtained from an extended UML 
fit to the two-dimensional AE-M\,c distribution. The 
projections of the fit results into AE and Afbc are 
shown in Fig. 7 and the fit results are summarized in 
Table II. The fit yields a signal of (136 ± 17) B+ 
(j)K+-/ and (35 ± 8) S° ^ 0X^7 candidates. The 
signal in the charged mode has a significance of 9.6 a, 
whereas that for the neutral mode is 5.4 tr, including 
systematic uncertainties. 

We also search for a possible contribution from 
kaonic resonances decaying to c/jK. To unfold the 
M^K distribution, we subtract all possible back- 
grounds and correct the (f>K invariant mass for the 
efficiency. The background-subtracted and efficiency- 
corrected M^K distributions are shown in Fig. 8. 
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Figure 7; The AE and Mbc projections for B"*" — > (f>K'^'y (upper) and B" — > 4>Kgj (lower). The points with error bars 
represent the data. The different curves show the total fit function (solid, red), total background function (long-dashed, 
black), continuum component (dotted, blue), the 6 — > c component (dashed-dotted, green) the non-resonant component 
as well as other charmless backgrounds (filled histogram, magenta). 



Table II The signal yields (Y), corrected efficiencies (e), 
branching fractions (B) and significances (<S) for the 
cpK^j and (jiK'^'y decay modes. 



Mode 


Y 


6(%) 


B (10-^) 


5(a) 




136 ± 17 


15.3 ±0.1 


2.34 ±0.29 ±0.23 


9.6 


(l>K°'y 


35 ±8 


10.0 ±0.1 


2.66 ±0.60 ±0.32 


5.4 



Nearly 72% of the signal events are concentrated in 
the low-mass region (1.5 < M^k < 2.0 GeV/c^). It is 
clear that the observed (jjK mass spectrum differs sig- 
nificantly from that expected in a three-body phase- 
space decay. 



7. Summary 

In summary, we report the first observation of ra- 
diative B'^ 'pKg'y decays in Belle using a data sam- 
ple of 772 X 10^ BB pairs. The observed signal yield 
is (35 ± 8) with a significance of 5.4 cr including sys- 
tematic uncertainties. We also precisely measure the 
B~^ <pK~^^ branching fraction with a significance 
of 9.6 tJ. The signal events are mostly concentrated 
at low (pK mass, which is similar to a two-body ra- 
diative decay. We report the measurements of time- 
dependent CP violation parameters in B^ K^tt'^j 
and B° K%p°-f decays, using 535 x 10^ and 657 x 10^ 
BB pairs, respectively. The Heavy Flavor Averaging 
Group (HFAG) [14] summary of the measurements of 




1.5 2 2.5 3 3.5 4 

U(i> K) (GeV/c^) 

Figure 8: The background-subtracted and efficiency- 
corrected <j)K mass distributions for _B+ — > (pK'^'y (upper) 
and B" — > (pKg'y (lower). The points with error bars rep- 
resent the data. The yield in each bin is obtained by the 
fitting procedure described in the text. The three-body 
phase-space model from the MC simulation is shown by 
the circles (blue). 



the parameter S in b sj modes by the Belle and 
BaBar is shown in Figure 9. With the present statis- 
tics, these measurements are consistent with the stan- 
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Figure 9: Summary of Belle and BaBar measurements of 
b ^ S"f Scp ■ 



dard model predictions and there is no indication of 
New Physics from right-handed currents in radiative 
B decays. The neutral B'^ 4'Kg'y mode has enough 
statistics for a future measurement of time-dependent 
CP violation. More luminosity is necessary for a pre- 
cise test of the SM. 
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